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Abstract - We study the nonequilibrium dynamics of multiband BCS superconductors subjected 
to ultrashort pump pulses. Using density-matrix theory, the time evolution of the Bogoliubov 
quasiparticle densities and the superconducting order parameters are computed as a function 
of pump pulse frequency, duration, and intensity. Focusing on two-band superconductors, we 
consider two different model systems. The first one, relevant for iron-based superconductors, 
describes two-band superconductors with a repulsive interband interaction V12 which is much 
larger than the intraband pairing terms. The second model, relevant for MgF32, deals with the 
opposite limit where the intraband interactions are dominant and the interband pair scattering V12 
is weak but attractive. For ultrashort pump pulses, both of these models exhibit a nonadiabatic 
behavior which is characterized by oscillations of the superconducting order parameters. We 
find that for nonvanishing V12, the superconducting gap on each band exhibits two oscillatory 
frequencies which are determined by the long-time asymptotic values of the gaps. The relative 
strength of these two frequency components depends sensitively on the magnitude of the interband 
interaction V12. 
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The discovery of iron-based superconductors [T] and 
MgB2 [2], where superconductivity is characterized by 
more than one order parameter, has lead to a renewed 
interest in multiband superconductors. The existence of 
multigap superconductivity was predicted more than fifty 
years ago in the pioneering works by Moskalenko [3] and 
Suhl et al. [3] , where it was argued that a band-dependent 
interaction potential can give rise to different order param- 
eters on different Fermi surface sheets. The presence of 
multiple superconducting order parameters has many in- 
teresting physical consequences, most of which are related 
to the relative amplitudes and phase differences between 
the gaps on different Fermi surfaces. 

In recent years, numerous studies of the nonequilib- 
rium dynamics of multiband superconductors have been 
performed using femtosecond time-resolved spectroscopy 
[5Ul0j. The relaxation kinetics measured in these experi- 
ments gives important information on the electronic band 
structure, on the electron-phonon coupling strengths, as 
well as on the symmetry of the superconducting order 
parameters. For example, time-resolved measurements 



on the iron pnictide superconductor Bao.6Ko.4Fe2As2 |10| 
have revealed two distinct quasiparticle relaxation pro- 
cesses: a fast one, whose decay rate increases linearly with 
excitation density, and a slow one, which is independent 
of excitation density. This behavior has been attributed 
to the multigap structure of the superconductor. Further- 
more, from a careful analysis of the temperature depen- 
dence of the decay rates, the authors of Ref. [TU] have 
been able to deduce the pairing symmetries of the super- 
conducting state. 

In this paper, we shall not be concerned with the relax- 
ation dynamics, but rather with the nonequilibrium evo- 
lution of multiband superconductors at times shorter than 
the relaxation time, i.e., with the nonadiabatic behavior of 
these system. It has been shown that the nonequilibrium 
response at these short time scales exhibits interesting 
properties, such as order parameter oscillations [TT1 - [T4] . In 
a pioneering work, B. Mansart et al. [T3] have been able to 
detect these coherent oscillations of the order parameter 
in a high-T c cuprate superconductor via transient broad- 
band reflectivity measurements. The recently discovered 
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multiband superconductors MgB2 and iron pnictide su- 
perconductors, which all have a relatively large supercon- 
ducting gap, offer an excellent opportunity to analyze this 
nonadiabatic dynamics also in multiband systems, where 
the interplay of the intraband pairing interactions and in- 
tcrband pair scattering terms play an important role. One 
goal of the present work, is to provide guidance for exper- 
imental studies of these interesting nonadiabatic response 
properties in multigap superconductors. 

In the folowing, we employ the density-matrix theory to 
study the temporal evolution of Bogoliubov quasiparticlc 
distributions in multiband superconductors after excita- 
tion by a short pump pulse. We focus on the nonadiabtic 
regime, i.e., the regime where the pump pulse duration t 
is much shorter than the dynamical time scales of the su- 
perconducting order parameters. In this case, the pump 
pulse creates rapid oscillations in the quasiparticle densi- 
ties. The essential physics of multiband superconductors 
can be captured by a minimal two-band model with Hamil- 
tonian H = H + H int , with H = J2 k i c £ M a tl<x a ida and 
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where a\^ a is the electron (hole) annihilation operator with 
band index I = 1,2, spin a, and momentum k. The bare 
dispersions of the two bands, £ki and £k2, are assumed 
to be quadratic with em = h 2 k 2 / (2m;) — E-pi, where m; 
and Epi represent the effective electron (hole) mass and 
Fermi energy of band I, respectively. The intraband pair- 
ing interactions are denoted by V££, and V££,, while the 
interband pair scattering is represented by V^£, = (V£} k )*. 
Assuming that the splitting of the bands is much larger 
than the superconducting energy scales, we neglect all in- 
terband pairing terms, i.e., we set (a_kzj,akZ't) = for all 
I 7^ I' . Furthermore, we restrict ourselves to momentum- 
independent pairing interactions with V££, = Vn, V££, = 
V22, and V^, = V12 = \V~i2\e 10 . The intraband pair cou- 
plings are assumed to be attractive, i.e., Vn, V22 > 0. The 
phase 9 of V12 is taken to be either or ir, depending 
on whether the interband pair scattering is attractive or 
repulsive. 

The dynamics of superconductors at times shorter than 
the quasiparticle energy relaxation time can be fully de- 
scribed within mean- field BCS theory [HJIT2]. We there- 
fore perform a mean-field decoupling of the quartic inter- 
actions in Eq. (|2|), which yields the BCS Hamiltonian 



#bcs = H a + ( A * a Lt a -ku + A ;* a ~kii a k 
kew,! 

with the mean-field gaps 

Ai = (Vn(a_ki4.fflkif) + ^i2(a-k24.ak2t)) , 
kew 



(2) 
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(V22(a-k2j.ak2t) + V2i(a_ki4.akit)) 



(3) 



The sums in Eqs. ([2]) and ((3]) are over the set W of mo- 
mentum vectors with \eyu\ < Tilu c and w c the cut-off fre- 
quency. The quasiparticle spectrum of Hbcs is given by 
{-E kl ,-Ete,+E kl ,+E k2 }, where E kl = ^4z + |A ; | 2 , 
with I = 1, 2. In the equilibrium, the order parameters Ai 
and A2 are determined by the self-consistency equation 
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where K t = £ keVV tanh (E M /2T) /E M , with 1 = 1,2. The 
gap equation ^ has a nontrivial solution whenever the 
determinant of the corresponding secular matrix vanishes, 



i.e., 
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From this condition, the transition temperature T c can be 
obtained. 

The two-band superconductor @ is perturbed by a 
short Gaussian pump pulse, which in the Coulomb gauge 
is described by a transverse vector potential A q (i) = 

iwot+x p+'"o') , with amplitude 
qo&x ! and full width 
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Aoe y , photon wave vector qo 
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at half maximum r. The optical pump pulse is centered 
at t = and has a central frequency u>o- The coupling of 
the vector potential A q to the two-band superconductor 
is given by 



H em (t) = eh £ „t 
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By means of the density-matrix formalism, we numer- 
ically compute the coherent response of the above model 
Hamiltonian after excitation by a short pump pulse. To 
this end, we need to derive equations of motion for the Bo- 
goliubov quasiparticle densities. It is advantageous to per- 
form this calculation in the basis in which -Hbcs 1 Eq. (|2|) , 
at t = ti, i.e., in the initial state, takes diagonal form. 
Hence, we consider the following Bogoliubov transforma- 
tion 



kz^kz 
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kZ> "-kZ| 
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where the Bogoliubov coefficients Ukz and «kz are time- 
independent, with Ukz = yl/2(l + Ski/E^i) and Vkz = 

e^VV2 0-Eki/Eu), where E kl = + I MU) \ 2 - 

The phase factors of the superconducting gaps in the ini- 
tial state = A;(t,)/ | Ai(ti) | are chosen to be either 
(+1)' or (—1)', depending on whether V12 in the initial 
state is attractive or repulsive. Substituting Eq. (J7J) into 
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Fig. 1: (Color online) Calculated time evolution of |A CT (£)| and 
|A„(i)| for the case of MgB 2 , for V12 = 0.2Vn and three differ- 
ent pump pulse widths r = 0.4 ps (solid green), 1.0 ps (dotted 
red), and 5.0 ps (dashed blue). The central energy of the pump 
pulse is huiQ — 8.0 meV, i.e., it is slightly larger than 2A 7r (t;), 

but smaller than 2A rT (t i ). We take A\t « 5.3 x 10~ 27 J s 



Eq. © yields 
#bcs(*) = 
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Note that here, as before, all terms except Aj(t) are time- 
independent. In particular, we have assumed that £ki has 
negligible time-dependence. Combining Eqs. ^ and © 
we find that the superconducting order parameters at time 
t can be expressed as 



Ai(t) 



E 

kew,; 



Wki'^aL'>(*)-«ki'<«ki'Ad'X*) 
Vw' ((4r«kr)W + (/3 k rV)W-l)]-(9) 



Hence, the temporal evolution of the gaps is fully de- 
termined by the time-dependence of the Bogoliubov 
quasiparticle densities (o^ot^, j) , (Ph^n), ( a kz<^k'z)' an< ^ 
(<^ki a k'i)- ^ s s h° wn m the Appendix, it is straightforward 
to derive a closed set of equations of motion for these quan- 
tities, see Eq. (fTTj) . 

We numerically solve the equations of motion for the 
Bogoliubov quasiparticle densities (see Appendix) to de- 
termine the time evolution of the order parameter ampli- 
tudes |Aj(f)|. As is evident from Eq. (JTTJ) , off-diagonal el- 
ements in the quasiparticle density matrices, such as, e.g., 
( a L a k+nq i' > ' are °f or der |Ao|™- Thus, for sufficiently 
small |Aof, the off-diagonal entries decrease rapidly as n 
increases. Therefore we set all off-diagonal elements with 



n > 4 to zero, which substantially reduces the computa- 
tional costs. For the numerical simulations we choose the 
BCS ground state at zero temperature as the initial state. 
In what follows, we consider two model systems, one with 
weak attractive V12, relevant for MgB2, and the other with 
strong repulsive V12, relevant for iron-based superconduc- 
tors. 

We start by discussing the model appropriate for the 
multiband superconductor MgB2. The electronic struc- 
ture of MgB2 consists of two sets of bands: cr-bands 
forming quasicylindrical hole pockets and 7r-bands giv- 
ing rise to tubular networks of electronlike Fermi surface 
sheets. Superconductivity in MgB2 arises from couplings 
between the electrons and the E2 9 -phonon mode. These 
couplings lead to attractive intraband interactions, which 
are stronger in the cr-bands than in the 7r-bands. Since 
the interband pair scattering is weak, this gives rise to two 
gaps with different magnitudes on the a- and 7r-bands. 

We describe the band structure of MgB2 in a somewhat 
oversimplified manner by two parabolic bands |16j . i.e, a 
hole band, with Fermi velocity Vp = 0.273 x 10 6 m/s and 
effective mass m a = — 3mo, and an electron band, with 
Vp = 1.0 x 10 6 m/s and = mo. Here, too denotes 
the free electron mass. The intraband pairing interac- 
tions on the a- and 7r-bands are chosen to be V\\ = 24.8 
meV and V22 = H-6 meV, respectively, whereas the at- 
tractive interband coupling is taken to be much smaller 
than either Vn or V22, i-C, we set V12 = 0.2Vn. With 
these material parameters, the superconducting gaps in 
the initial state Aj(tj) can be computed from the self- 
consistency equation (j4]). Assuming that the cut-off en- 
ergy huj c is approximately equal to the E2 9 -phonon energy, 
i.e, huj c = 50 meV, we find that in the initial state the gaps 
on the a- and 7r-bands are given by A CT (t.;) = 7 meV and 
A 7r (ti) = 3 meV, respectively. The central energy of the 
optical pump pulse hujQ is chosen to be either in between 
2A 7r (fi) and 2A CT (i.;), or equal to twice the gap on the 
er-band, i.e., Tiujq = 14 meV. 

Let us first discuss the most interesting case, namely 
the situation where the pump photon energy fiojQ lies in 
between 2A 7r (ii) and 2A CT (^). In Figs. [T^a) and[ljb) we 
plot the calculated time evolution of |A CT (t)| and |A„.(f)|, 
respectively, for huj = 8 meV, V12 = 0.2Vii, and three 
different pump pulse lenghts t = 0.4, 1.0, and 5.0 ps. 
In accordance with previous results on single band super- 
conductors [T2"Hl~4"l [T7] . we observe two distinct dynam- 
ical regimes: a nonadiabatic regime for pulse widths r 
much smaller than the dynamical time-scales of the su- 
perconductor r A ~ h/(2\Ai\) and an adiabatic regime 
for r 3> t A; . In the nonadiabatic regime the Bogoliubov 
quasiparticle densities build up coherently while the sys- 
tem is out of equilibrium. This gives rise first to a mono- 
tonic increase and then to fast oscillations in the quasi- 
particle occupations. Similarly, the superconducting gaps 
first decrease monotonically and then start to oscillate 
rapidly (solid green curves in Fig.[T]). Since our model sys- 
tem does not contain any relaxation terms, such as inter- 
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actions or couplings to an external bath, the oscillations in 
the quasiparticle occupations are undamped. In contrast, 
the order parameter oscillations decay as 1/%/i. That is, 
due to decoherence effects, the superconducting gap Aj(t) 
approaches a constant value as time tends to infinity. To 
be more precise, we find that to a good approximation the 
time evolution of |Aj(i)| is described by 

i 2 

|A,(t)| = Ar + b n' cos ( 2A ™t + <t>w) , (10) 

where bw and <pui are constants that depend on the ini- 
tial conditions. In other words, the superconducting gaps 
oscillate with an amplitude decaying as 1/ \ft and two fre- 
quencies, which are determined by the long-time asymp- 
totic gap values Af° < |A;(ij)| [see Figs. Ufa) andQJb)]. 
Interestingly, the relative strength of the two frequency 
components wa„ = 2A^°/fi and oja x = 2A^/h depends 
sensitively on the magnitude of Vi2- That is, the ampli- 
tudes &12 and 621 increase with increasing intcrband pair 
scattering V12. For V12 = 0, on the other hand, the am- 
plitudes 612 and 621 vanish, and hence Eq. (fTU| reduces to 
the well-known result for single band models, where the 
gap oscillates with a single frequency P^HT^llTTlfTB] . 

A nonvanishing interband pair scattering V12 not only 
affects the relative strength of the oscillation frequencies 
but also strongly influences the long-time asymptotic gap 
values A^° and hence the frequencies of the order param- 
eter oscillations. This is illustrated in Figs.^a) and[5Jb), 
which show the dependence of A; 00 on the integrated pump 
pulse intensity A\r for Tiwo = 8 meV and both for zero 
and nonzero V12. Comparing Figs. [3Ja) and^b), we ob- 
serve that in the presence of nonzero interband pair scat- 
tering the time evolutions of the two gaps are coupled 
together. The change in A z °° due to interband interac- 
tions seems to be roughly proportional to V^Af?, with 
I' ^ I [cf. Eq. ©]. That is, the larger gap A a {t) influ- 
ences the smaller gap A„ (t) more strongly than vice versa. 
With increasing integrated laser intensity both gaps are 
suppressed simultaneously, until for sufficiently large A\t 
they vanish at the same value of A\t. For zero V\%, in 
contrast, the order parameters can be suppressed inde- 
pendently. For example, we find that in the nonadiabatic 
regime the larger gap A CT (t) vanishes at a smaller inte- 
grated intensity than the smaller gap A„(t) [solid green 
curve in Fig.[2Jb)]. We observe that in general the asymp- 
totic gap values Af° decrease linearly at small A%t, but 
deviate from this behavior at larger integrated intensities. 
With increasing AqT, the curves corresponding to short 
pump pulses (r = 0.4 ps) show a downward bend, while 
those with longer pump pulses (r = 1.0 and 5.0 ps) flatten 
due to Pauli blocking. 

Let us also briefly discuss the case where the pump pulse 
energy huo is equal to 2A (T (i i ) = 14 meV. In general, 
the behavior of the gaps in this case is rather similar to 
the previously discussed case, where Hujo = 8 meV. As 
before, we find that for pump pulse widths r 3> t Aj both 
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Fig. 2: (Color online) Asymptotic values of the order pa- 
rameters and as a function of integrated pump pulse 
intensity for three different pump pulse widths r. The central 
energy of the pump pulse is taken to be Tlljo = 8 meV in panels 
(a) and (b) , and huo = 14 meV in panels (c) and (d) . The left 
and right panels show the results for V12 = 0.2Vn and V12 = 0, 
respectively. The symbols represent calculated values, while 
the solid and dashed lines are a guide to the eye. Here, the 
arbitrary units in Aqt have to be multiplied by the constant 

-97 1 2 S 3 

d w 1.06 x 10 ^ 3 to get the physical units. 

gap amplitudes | A; (t) | decrease monotonically from their 
initial equilibrium values to their final values A^° , whereas 
for t <C t A( the gaps approach their asymptotic values in 
an oscillatory fashion according to Eq. (fTU|) . However, 
as can be seen from Figs. [2jc) andjSJd), the dependence 
of Aj°° on the integrated intensity AqT is quite different 
than in Figs. [5Ja) and[^b). Interestingly, we find that a 
pump pulse with central energy Tiujq — 2A a (ti) depletes 
the superconducting condensate more efficiently on the a- 
band than on the 7r-band. In particular, for V12 = the 
gap on the 7r-band is almost unaffected by the pump pulse, 
even for large A^t, whereas the asymptotic gap value on 
the CT-band decreases steadily with increasing A\t. For 
nonzero V12 the behavior is similar, although here the gap 
on the 7r-band is slightly suppressed, which is due to the 
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Fig. 3: (Color online) (a) Calculated time evolution of |Ah(t)| 
and |A e (i)| for the case of iron-based superconductors. The 
inset shows the spectral distribution of the gap oscillations for 
the same parameters as in the main panel, (b) Asymptotic 
values of the order parameters A^° and A£° versus integrated 
pump pulse intensity for three different pump pulse widths r. 
In this figure we choose V12 = — 10 Vn and Tiujq = 14 meV. In 
panel (a), we set r = 0.4 ps and A\t w 1.44 x 10~ 26 Js/Cm 
[corresponding to the point Aqt = 13 arb. units in panel (b)]. 
The symbols in panel (b) represent calculated values, while the 
lines are a guide to the eye. Here, the arbitrary units in A\t 

— 27 T 2 S 3 

have to be multiplied by the constant d w 1.11 x 10 c' 2 2 
to get the physical units. 



coupling with the gap on the er-band. 

Secondly, we consider a model appropriate for the iron 
pnictide superconductors [T5]. Many features of iron- 
based superconductors can be captured within a simple 
two-band model with quasinested holelike and electron- 
like pockets centered at the T and M points in the Bril- 
louin zone, respectively [20]. Due to the quasinesting of 
the Fermi surfaces, the repulsive interband pair scatter- 
ing in such a description is much larger than the intra- 
band pairing potentials. The latter are assumed to be 
overscreened due to the nesting effects and become at- 
tractive for low energies relevant for superconductivity. 
Therefore, in what follows, we model the band structure 
of iron-based superconductors by two parabolic bands, a 
hole band with effective mass = — 3toq and Fermi 
velocity Up = 0.168 x 10 6 m/s, and an electron band 
with m c = m and up = 1.453 x 10 6 m/s [2T]. For 
convenience we shift the location of the electron Fermi 
surface from the M point to the T point, i.e., we set 
£2k+Q = £2k, where Q = (tt, 7t). The attractive intraband 
pairing interactions on the electron and hole Fermi sur- 
faces are chosen to be Vn = 2.466 meV and V22 = 2.463 
meV, while the repulsive interband pair scattering is set to 
V12 = —24.42 meV. With these material parameters and 
assuming that the cut-off energy hui c is determined by the 
energy of the spin fluctuations in the paramagnetic state, 
hio c = 30 meV |22j . we find that the equilibrium gaps are 
given by Ah(ti) = —A e (ti) = 6.5 meV. That is, the gaps 



on the hole and electron bands are equal in magnitude but 
opposite in phase, which is roughly in line with experimen- 
tal findings in some of the hole-doped iron-based supercon- 
ductors |21j . The two-band superconductor is driven out 
of equilibrium by a pump pulse with energy Hluq = 14 
meV, which is of the same order but slightly larger than 
twice the gap amplitudes in the initial state. 

Based on this model for pnictide superconductors, we 
compute the temporal evolution of the superconducting 
order parameters on the hole and electron bands. In 
Fig. G3a), the time dependence of |Aj(t)| is shown for a 
pump pulse with length r = 0.4 ps, i.e., for the nona- 
diabatic regime. We observe that both gaps exhibit al- 
most the same behavior: they first decrease monotoni- 
cally and then start to oscillate with different frequen- 
cies producing a beatinglike pattern in the amplitudes. 
Due to the smaller Fermi velocity and the larger effec- 
tive mass of the hole band, the gap on the hole band 
Ah(i) is slightly less suppressed than the gap on the 
electron band A c (i). Because of the nonzero interband 
coupling V12, both gaps oscillate with two frequencies, 
uj Ac = 2A™/h and w Ah = 2Ag°/7i, which differ very lit- 
tle, i.e., \oJA h — W A C | "C wa c . This in turn leads to a pro- 
nounced beating phenomenon, as can be seen in Fig.[3ta). 

The dependence of the asymptotic gap values A^° on 
the integrated pump pulse intensity A\t is plotted in 
Fig. [3jb) . Due to the large interband coupling V12 the 
gaps on both bands decrease almost synchronously with 
increasing AqT. However, the asymptotic value of the gap 
on the electron band A£° is always slightly smaller than 
the asymptotic value of the gap on the hole band A£°. 
As discussed above, this is because of the larger effective 
mass and the smaller Fermi velocity on the hole band. In 
general, we find that the relative difference between A£° 
and AJ?° gradually increases with increasing pump pulse 
intensity. Correspondingly, the difference between the two 
oscillation frequencies WA h and WA h increases, and hence 
the beating phenomenon becomes less pronounced, as the 
integrated pump pulse intensity is increased. 

In this work we analyzed the nonequilibrium dynamics 
of two-band superconductors after excitation by a short 
opitcal pump pulse. We considered two model systems: 
one with dominant intraband pairing and weak attractive 
interband pair scattering, and one where the repuslive in- 
trerband interactions are much larger than the intraband 
pairing potentials. The former model is relevant for MgB2 
superconductors, whereas the latter one is appropriate for 
iron-based superconductors. For both of these model sys- 
tems we numerically computed the time evolution of the 
order parameters as a function of pump pulse duration and 
integrated pump pulse intensity. Our main observation is 
that the ratio between the gaps for asymptotically large 
times depends sensitively on the interband Cooper-pairing 
strength and differs from its equilibrium value. This al- 
lows in principle to use pump-probe experiments for di- 
rect identification of the interband Cooper-pair scattering 
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strength. In addition, sufficiently short pump pulses cre- 
ate fast oscillations in the gap amplitudes and the quasi- 
particle densities of these two-band superconductors. We 
found that for nonzero interband pair scattering V12 these 
oscillations are characterized by two frequencies which are 
determined by the long-time asymptotic values of the gaps 
on the two bands (see Fig. [TJ . We showed that the relative 
strength of the two frequency components sensitively de- 
pends on the magnitude of the interband interaction V12 ■ 
When the gaps on the two bands are of similar magnitude 
(which is the case, for example, in iron-based supercon- 
ductors) the relative difference between the two oscillation 
frequencies is small, and hence the gaps oscillate with a 
beatinglike pattern in the amplitudes. 
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Appendix. In this Appendix, we present one of the 
equations of motion (as an example) which determine 
the time evolution of the quasiparticle densities (o^ct^,,), 

(/^k/^k'/)' ( a kiAc';)> and (pii a t>i)- Here > as m tiic main 

text, l~l,2 denotes the band index. Since (a k ;/3 k /;) and 
(^k'i a Ll) are related by Hermitian conjugation, only six 
out of the eight quasiparticle densities need to be evalu- 
ated numerically. By use of Heisenberg's equation of mo- 
tion, the following equation for (a k ;a k , ; ) can be obtained 

-?/2k'i(/?iUaL) + ?72kz("k';Ad) - ^ k ■ A i 

q=±qo 

( M *k',k'-q( a L"k<-q() + ^k'.k'-q^k'-qjaL) 
-^k+q,k( Q k+qi a k'i) - ^ktk+q^k'iAc+q*)) 

+ ^E( E Aq_q,-A q/ ) 

q \q/=±qo / 

( M ik',k<-qKz a k'-qi) + ££'_ q ,k'<£k'-ql Qt k|) 
- M ik+q,k(«k+qi«k'i> " £i* k+q (ak'zAc+qi)), 

(11) 

where we have introduced the short-hand notation 

M ztk' = U ki"k'Z ± V W<>1 and L Zk,k' = W kZ«k'Z ± WkZUk'Z, 

with the coherence factors Ukz and t>kz. The functions 
??ikz and ?72kz are defined in the main text. In deriving 
Eqs. (fill , we used the relation A* = A_ q . 
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